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ADAPTATION OF FLUX-CORRECTED TRANSPORT ALGORITHMS 
FOR MODELLING BLAST WAVES 


Blast wave phenomena include reactive and two-phase flows associated 
with the motion of chemical explosion products; the propagation of shocks, 
rarefaction waves, and contact discontinuities through a nonideal medium 
(real air, possibly thermally stratified and containing dust and water 
vapor); and the interaction of the blast waves (including boundary layer 
effects) with structural surfaces. Flux-Corrected Transport (FCT) represents 
an accurate and flexible class of methods for solving such nonsteady 
compressible flow problems (Boris and Book, 1976 ). Coupled with a 
nondiffusive adaptive gridding scheme (Book, et al ., I98O; Fiy, et al ., 

1981), it enables complex time-dependent shocks to be efficiently 
"captured." 

In models which treat all the physical effects required for blast wave 
simulation, truncation errors inherent in the underlying finite-difference 
scheme are exacerbated by nonlinear coupling between the fluid equations and 
by the greater conylexity of the phenomena being simulated. Typical of these 
errors are the "terraces" which develop under some circumstances on the 
flanks of sloping profiles when the growth of ripples due to phase errors at 
short length scales is terminated by the action of the flux limiter. Two 
approaches are possible toward eliminating them: inproving the short- 
wavelength phase and amplitude properties of the underlying algorithm, and 
switching on additional diffusion locally. The latter approach folds 
information about the shape of the profile and the nature of the physical 
process taking place (e.g., rarefaction) into the switch criterion, thus. 
changing the FCT technique from a "convective equation solver" to a "fluid 
system solver." In doing this, care nust be taken to avoid losing the 
accuracy, robustness and problem-independence which constitute valuable 
attributes of FCT algorithms (Book, et al ., 1981). 

Tests carried out on scalar advection of simple density profiles by a 
uniform flow field show that terracing does not require either diverging 
velocities or discontinuities in the profile, but appears topically (for 
0) where the first and second derivatives of density have the same sign 
(Fig.l). In order to improve the properties of the basic difference scheme, 
we propose a new algorithm for integrating generalized continuity equations 
over a timestep < 5 t. Consider the following three-point transport scheme. 
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The arrays {p?} and {p 1 ?} are the old and new densities, p. and p. are 

temporary intermediate densities, and n, 6, k, X, and y are velocity- 
dependent coefficients. Here k and X are diffusion coefficients, and y is 
the antidiffusion coefficient. In the actual algorithm, is corrected 

1 .r* . n_ 


(hence the name FCT) to a value <)£ 


J +1 /g 


chosen so no extrema in pj can be 


enhanced or new ones introduced in pj. Previous FCT algorithms had 9=0; 

the widely used ETBFCT and related algorithms (Boris, 1976) have in addition 

< = 0. If we define p^ to be sinusoidal with wave number k on a mesh with 

uniform spacing Sx, so that p°? = exp (ij 3) where 3 = k6x, then the new 

J 

density array satisfies 

p^/p? = A = l-2i(n+9)sin3 +2( k+X) (cos 3-1) 

J J 

- 2y(cos 3-1) [l-2insin3+2ic(cos 3-l) ]. 


From A we can determine the amplification a = A and relative phase error R = 
(1/e3)tan - 1 (-ImA/ReA)-l, where e = v 6 t/ 6 x is the Courant number. 

Expanding in powers of 3 cue find 

a = 1 + ag3^ + 0^3^ + + • • • ; 

R = Rq + R23^ + Rj^3^ + Rg3^ + . . . . 

First-order accuracy entails making Rq vanish, which requires that n + 0 = 
e/2. Second-order accuracy (oig = 0) implies that y = k + X - e 2 /2. 
Analogously, the "reduced-phase-error" property Rg = 0 (Boris and Book, 

1976) determines y = (l-e 2 )/ 6 , thus leaving two free parameters. One of 
these can be used to make R^ vanish also. The resulting phase error R(3) 
is small not only as 3^0, but also for larger values of 3 , corresponding to 
the short wavelengths responsible for terraces (Fig. 2). The remaining 
parameter n can be chosen to relax the Courant number restriction needed to 
ensure positivity from e < 1/2 to e < 1. When coded, these changes 
necessitate a small increase in the operation count of ETBFCT along with a 
small increase in overhead to precalculate the two new arrays of velocity- 
dependent transport coefficients. On advection tests, the new algorithm 
completely eliminated terraces (Fig. 3). When applied to the coupled systems 
of gas dynamic equations, it produced profiles which closely approximate the 
Riemann solution of the exploding diaphragm problem (Fig. 4). 

The second approach uses a rarefaction flux limiter (RFL) to eliminate 
numerical ripples in strong rarefaction waves. This approach is physically 
motivated. Raw anti-diffusive fluxes <f ) ^ + T /2 are limited so that 
the slope of local flow field profiles decays with time in a rarefaction 
wave. In effect, additional diffusion is left in the field to maintain 
monotonicity of local slopes. For multi-material calculations a "contact 
surface sensor" is needed to detect physical discontinuities and shut off the 
RFL locally. 

In addition we found that some care was required when applying 
generalized continuity equation solvers to a system of equations. Truncation 
errors of the various equations can interact, causing undershoots or 
overshoots in nonconvective quantities such as pressure. We found that it 
was necessary to monotonize derived quantities (pressure, velocity) before 
using them in minimal-diffusion transport algorithms. 
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The above methodology has been applied to a series of test problems 
initiated by a spherical high-explosive (HE) detonation in air. An ideal 
Chapman-Jouguet detonation was used to specify the initial conditions; 
afterburning was neglected. In the absence of reflecting surfaces, 
spherical symmetry is maintained and the calculation remains one¬ 
dimensional. A nonuniform radial grid was used with extremely fine zoning 
near the shock front. The grid was moved so that the shock remained 
approximately fixed with respect to the mesh. The original version of the 
FCT algorithm gave rise to pronounced terraces in the rarefaction region. 

This would have rendered any two-dimensional calculation involving shock 
diffraction or nonideal effects dubious. The techniques described here 
improved the blast wave results considerably. The decrease in phase error 
reduced terracing dramatically. 

Next, a two-dimensional (2D) numerical calculation was performed to 
simulate one of Carpenter's (197*0 height-of-burst experiments which used 
spherical 8-lb. charges of PBX 9404 at 51*6 cm. The previous fine-zoned ID 
calculation was used to initialize the problem. It was mapped onto the 2D 
grid just prior to the onset of reflection. The solution was then advanced 
in time, with pressure being calculated from a real-air equation of state and 
a JWL equation of state for the combustion products. The front of the blast 
wave was captured in a finely gridded region which moved outward horizon¬ 
tally. Special care was taken to ensure that the grid moved smoothly. The 
resulting solution, particularly the curve of peak overpressure vs. range, 
was consistent with Carpenter's experimental data (Fig. 8). Although this 
calculation represents a reasonable accurate simulation of the double-Mach- 
stem region, no doubt improvements can and will be made to numerically model 
such phenomena. 
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0.0 RADIUS - CM 

(a) 



(b) 

Fig. 1 — Rounded half circle used in passive scalar advection tests (a) initially, 
and (b) after propagation for 14 cycles using JPBFCT. Note that terraces form 
even, as here, in the absence of corners in the profile. Tick marks indicate 
computational zones (N = 100). 
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HULTICDEFFICIENT 


REL. PHASE 



Fig. 2 — Contour plot of R(j3,e) for new multicoefficient FCT algorithm. 
Note R » 0 except for 0 > 3 tt/ 2. The relative phase error vanishes ex¬ 
actly for e = 1/2 and e = 1. 
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JPBFCT ROUNDED HALF-CIRCLE 
CYCLE = 14 TIME = 2.17 10 3 



(a) 


LCPFCT ROUNDED HALF-CIRCLE 
CYCLE = 14 TIME = 2.17 10 



(b) 

Fig. 3 — (a) Blowup of Fig. 1(a) (dashed line) compared with (b) same profile as 
computed using new sixth order-phase-accurate FCT algorithm. Solid traces are 
exact solutions. 
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(b) 

Fig. 4 — (a) Exact and (b) computed solution of exploding diaphragm problem 
(10-to-l initial density jump, 100-to-l initial pressure jump) 
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Fig. 5 — One-dimensional solution of expanding HE products and air calculated with 
the new algorithm using 500 equally spaced zones. Note contact surface separating 
HE products from air. 
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ZONES 60 COARSE ZONES 




Fig. 6 — Adaptive grid for height-of-burst problem shown (a) initially and 
(b) at time when transition to Mach reflection occurs 




















PMAX - 1.74*11 1 TIME- S.21*10' 

PRESSURE VS TIME 

Fig. 7 — Pressure-time histories directly beneath burst site. Note second peak, associated with 
interaction between shock reflected from ground and following contact surface. 
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ire vs distance along ground surface. 
; Carpenter’s (1974) data. 
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(DESIRES NC CN'WDI) 
01CY ATTN AUL-LSE 


ASSISTANT CHIEF OF STAFF 

INTELLIGENCE 

DEPARTMENT OF THF AIP FORCE 
WASHINGTON* DC 201 PC 
01CY ATTN IN PM 4AR32 

ASSISTANT CHIEF CF STAFF 
STUDIES £ ANALYSES 
pfro^oTMENT OF THE AIR FO R CE 
WASHINGTON, DC 20320 

01CY ATTN AF/SAMI (TECH 119) 

ASSISTANT SECRETARY OF THF AF 
RESEARCH, DEVELOPMENT £ LOGISTICS 
DEPARTMENT OF THF AIR FCRCF 
WASHINGTON* DC 2033C 

01CY ATTN SAFALR/CEP FOR STRAT £ SPACE SYS 

BALLISTIC MISSILE OFF ICE/ON 
AIP FOPCE SYSTEMS CCMMAKT 
NO PTC* 1 AFS, CA R?4CR 
{VIN U T F MAN) 

01CY ATTN M NNXH G KALANSKY 

QICY ATTN VNMXH M CFLVECCHIO 

01CY ATTN V NN W CFABTPEE 

OICY ATTN YM m XH D GAGE 

0 1CY ATTN MMNX 
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0E°UTY CHIEF OF $TA P F 
RESEARCH* DEVELOPMENT, £ ACC 
DEPARTMEN t OF t F'F A IP FCFCF 
WASHINGTON, DC ?0?30 

OICY ATTN AFROOI N ALFXAND p OW 
OICY ATTN A r R D D N 
OICY ATTN AFPDQt 

DEoIJTY CHIEF OF STAFF 
LOGISTICS £ ENGINEERING 
nppjjpjMENT OF THE AI p FORCE 
WASHINGTON, DC 20 33C 
OICY ATTN LEFF 

COMMANDER 

FOREIGN TECHNOLOGY DIVISION, AFSC 
WRIGH T-oATTEP SON AFP, OH 45433 
OICY ATTN Ml IS L IPSARY 

COMMANDER 

poMf AIR DEVELOPMENT CENTER, AFSC 
GRTFFISS AF3, NY 134A I 

(C ESI n FS NC CNVDT) 

OICY ATTN TSLO 

STRATEGIC AIP COMMAND 
DEPARTMENT p F THF AI D FORCE 
OFFUTT AFB, MB 6«113 

OICY ATTN NPI-ST INFO LIBRARY 

OICY ATTN XPFS 

OICY ATTN INT J MCKINNEY 

VELA SEIS M OLOGICAL CENTEP 
312 MONTGOMERY STREET 
ALEXANDRA, VA 22 314 
OICY A T TM G OIL 0 ICH 

DE b APTMENT OF ENERGY 
AL BUG HER CUE O^CRATICAS OFFIC.F 
P 0 BOX 5400 
ALBUQUERQUE, N M P7115 
OICY ATTN CTID 

DEPARTMENT n F ENFPGY 
WASHINGTON, DC 20545 
OICY ATTN G’-LA/RDET 

DEPARTMENT OF ENERGY 
NEVADA OPERATIONS OFFICE 
P 0 BOX 14100 
LAS VEGAS, NV 30114 

OICY ATTN NAIL 6 PECOROS FOP TECHNICAL 


t I BPARY 
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LAWPENCF LIVE°MOPE NATIONAL LAB 
P 0 BOX BOP 
LIVERMORE, CA 94550 


OICY 

ATTN 

L-°0 R DCNG 



OICY 

V T T N 

L-205 J FEAPST 

(CLASS 

L-203) 

OICY 

ATTN 

L -90 o NORRIS 

(CLASS 

L-504) 

OICY 

ATTN 

L~? J KAHN 



OICY 

ATTN 

D GLENN 



OICY 

ATTN 

L 437 R SCHCCK 



OICY 

ATTN 

TECHNICAL INFO 

O p ?T . 

LIBRARY 

OICY 

ATTN 

L-200 T BITKOVICH 



L n S ALAMOS NATIONAL SCIFNTI p TC LAB 
V A I L ST A T I ON r -000 
P n BOX 1463 
LOS ALAMOS, MM P754 5 

(CLASSIFIED ONLY TC MAIL STATION 50C0) 
01CY ATTN R WH TT AK FR 
01CY ATTN C KELLER 
01CY ATTN MT SANFORD 
01CY ATTN 364 (CLASS REPORTS LIB) 

01CY ATTN E JONES 

LOVELACE BIOVFDICAL f. 

ENVIRONMENTAL PSCH INSTITUTE, INC. 

P 0 BOX 5890 
ALBUQ'JFRQIJF, N« B7I15 

OICY ATTN o JQN C S (UNCL CNLY) 

OAK PIDCE NATIONAL LABORATORY 
NUCLEAR DIVISION 
X-J.0 LAB RECORDS DIVISION 
P 0 BOX X 

CAK RIDGE, TN 37330 

OICY ATTN CIVIL DEF RES P90J 
OICY ATTN CENTRAL P.SCK LIBRARY 


SAND I A LABORATORIES 
LIVERMORE LABORATORY 
o o BOX 969 
LIVERMORE, CA 94550 

OICY ATTN LIBRARY € S C CU°IT' / CLAS SI r IC AT I CN DIV. 
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SV’DIA NATIONAL LAB 
P 0 POX 5800 
ALBUQUERQUE, N M S7195 

(ALL CLASS ATTN S C C CONTPOL OFC FOR) 


OICY 

ATTN 

A CHAD AN 

OICY 

A TT N 

l HILL 

OICY 

ATTN 

ORG 1250 W PRC 

OICY 

ATTN 

A CFA3 TA 

OICY 

ATT N 

W PfHERTY 

OICY 

ATTN 

3141 

OICY 

ATTN 

L VOPTMAW 

OICY 

ATTN 

J BENISTER 

MTP AL 

INTELLIGENCE AGENCY 


WASHINGTON» OC 205C5 
OICY ATT\i OSWP/NEC 


DEPARTMENT OF THF INTERIOR 
BUREAU CF MINES 

BLDG 20 « DENVER FFDEPAL CENTFP 
DENVER. CO P0225 

((UNCL CNLYJ) 

OICY ATTN TECH LIP (UNCI CNLY1 


DIRECTOR 

FEDERAL EMERGENCY MANAGEMENT AGENCY 
NATIONAL SEC OFC MITIGATION F. RSCH 
1725 I STREFT. NW 
V'ASHING T ON. DC 20A72 

((ALL CLASS ATTN S1C5 DCC CCNTRCL FCR)) 
OICY ATTN MTTTGATICN F RSCH DIV 

AC ' 1°E X C0°P• 

A«5 CLYDF AVFNM.JC 
MOUNTAIN VIEW, CA ° ADA ? 

OICY ATTN C WOLF 

AEPOS RACE COPP. 

P 0 BOX 92957 
LOS ANGELAS, CA 90009 
OICY ATTN H mirelS 

OICY ATTN TECHNICAL INFORMATION SERVICES 

ASBABIAN ASSOCIATES 
250 N NASH STREET 
El. SFGUNDO» C* 90?A5 
OICY ATTN * AGBA 8 IAN 

ANALYTIC SERVICES, TNC. 

AGO ARMY-NAVY OR TVE 
ARLINGTON, VA 22202 

OICY ATTN G HESSFLEACHER 
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APPLIED RESEARCH ASSOCIATES, INC 
?601 WYCMING BLVO NE SUITE F~1 
ALBUQUERQUE, f ' 1 ‘' B711? 

01CY ATTN J F R A T T 0 H 
OICY ATTN N FIGGING 

ADTLTED THED-Y, INC- 
1010 WESTWOOD PL VC 
LOS ANGFLES, CA 90"24 

(? CVS IF L NT LASS CP 1 CY IF CLASS! 
OIC.Y ATTN J TRULIC 

A° T EC ASSOCIATES, INC. 

26046 EDEN LANCING POA0 
HAYWARD* CA 94545 
OICY ATTN S GILL 

AGTRCN P5SAPCH £ ENGINEERING 
1901 OLD MIDDLFFIELD WAY 415 
•VC J N'T AI N VIEW, CA 5404 3 
OICY ATTN J HUNTINGTON 

A-/CD RESEARCH £ SYSTFVS CROUP 
201 LOWELI STREE T 
WIL M INGTON, MA 01837 

OICY ATT M LIBRARY AF30 

BDM CPPP. 

7915 JONES BRANCH DRIVE 
VCLFAN, VA 2 71 02 

OICY ATTN A LAVAC-KINC 
OICY ATJN T NEIGHPfPS 
OICY ATTN CORPORATE LIBRARY 

rhm rnpo, 
a n box 9274 
ALBUQUERQUE, \' M S7119 

OICY ATTN R HENSLEY 

DOE I ►'G CO. 

P 0 BOX 3707 
SEATTLE, WA 9R1?4 
OICY ATTN S STRACK 
OICY ATTN AFROS D ACE LIBRARY 
OICY ATTN M/S 42/37 P CARLSON 

CALIFORNIA 'RESEARCH £ TECHNOLOGY, INC. 

62 65 VAR TEL. AVENUE 
WOODLAND HILLS, CA 91367 
OICY ATTN LIBRARY 
OICY AT T u K KREYENHAGEN 
OICY ATTN M ROSENBLATT 



CALIFORNIA RESEARCH Z TECHNOLOGY, INC. 
4049 FIRST STREET 
LIVEPNORF, CA c 4550 
0ICY ATTN D PR^HAL 

CAL S»AN COR P . 
o 0 BOX 400 
BUFFALO, NY 14275 
OIC.Y ATTN LIBRARY 


OENVF D , UNIVERSITY rF 

COLORADO S E 9 IN A p V 

QFNVFR RFSFAPCH INST ITU T r 
P 0 POX 10177 
DENVER, CO 30210 

(ONLY 1 COPY C p CLASS R p T S ) 
01CY ATTN SEC OFFICER FC p J WISOTSKI 

EGCG WASH. ANAIVTICAL SVCS CTP, INC. 

D 0 FOX 1021P 
ALB'JQUEROUE» NM 87 114 
01CY ATTN LIBRARY 

FRTC H. WANG 

CIVIL ENGINEERING p SCH FAC 
UNIVERSITY OF NEW VFXICC 
UNIVFRSTTY STATION 
P 0 BOX 25 

ALBUQUERQUE, NY 87121 
01CY f ATTN J LAYS 
01CY ATTN P LCODE 
01CY ATTN N BAUM 
0 ICY ATTN J KGVARNA 

GARD, INC. 

744° N NATCHEZ AVENUE 
NIL FS, IL 60648 

0 ICY ATTN G N E10 H A T 0 T (UNCL ONLY) 

GENERAL ELECTRIC CO. 

S°ACE DIVISION 
VALLEY FQP.GE SPACE CENTER 
D 0 enx 8555 
PHI LADELPH I A, PA 19101 
0 ICY ATTN M 3QPTNEP 

GENERAL RESEARCH COPP. 

SANTA BARBARA DIVISICN 
P 0 BOX 6770 
SANTA BARBARA, CA 92111 
01CY ATTN TIO 


25 



H—TpCH LABS j INC. 

P 0 BOX 1686 
SANTA MONICA, CA OQ406 
01CV ATTN 8 HAP IE NT AIM 

HORIZONS TECHNOLOGY, INC* 

7830 CLAIPEMQNT MESA BLVD 
SAN 0 I EGO, CA °2111 
OICY ATTN R KRUGER 

11 T RESEARCH INSTITUTE 
10 W 35TH STREET 
CHICAGO, TL 60616 
01C.Y ATTN R WELCH 
01CY ATTN ,M JOHNSON 
01CY ATTN DOCUMENTS LIBRARY 

INPCPMAT ION SCIENCE, INC. 

123 W PAORF ST D EFT 
SANTA BARBARA, CA 93105 
01 C.Y ATTN W 01107 IAK 

INSTITUTE EOF DEFENSE ANALYSES 
400 APMY-NAVY DRIVE 
ARLINGTON, VA 22202 

01CY ATTN CLASSIFIED LIBRARY 

J D HALTIWANGER CONSULT ENG SVCS 
PM 106A CIVIL ENGINEERING BI.CG 
208 N ROM INF STREET 
"RSANAf IL 61801 
01CY ATTN W HALL 

J. H. WIGGINS CO., INC. 

1650 S PACIFIC COAST HIGHWAY 
REDONDO BEACH, CA 9C277 
01C Y ATTN J COLLINS 

KAMAN AVI DYNE 
83 SECOND AVENUE 
NORTHWEST INDUSTRIAL PARK 
BURLINGTON, MA 01^03 
0 1CY ATTN P PLETENIK 
01CY ATTN LIBRARY 
0 1CY ATTN N HOBBS 
0 1CY ATTN E CRISCICNE 

KAMAN SCI FNCFS COP''. 

° 0 BOX 7463 

COLORADO SPRINGS, CC 80933 
01CY ATTN D SACHS 
OICY ATTN F SHELTCN 
OICY ATTN LIBRARY 
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KAMAN TEMPO 

916 STATE S TP EE T ( D 0 DRAWER CQ) 

SANTA BARBARA, CA 92102 
01CY ATTN DASIAC 

LOCKHEED MISSILES S SPACE CC., INC. 

P 0 BOX 504 
SUNNYVALE, CA 94086 
01CY ATTN J WETSNER 
0 ICY ATTN TTC-LI BPARY 

MARTIN MARIETTA COR°. 

P 0 BOX 5937 
ORLANDO, EL 32P55 
OICY ATTN G ECTIFC 

MARTIN MARIETTA COPp. 
o 0 BOX 179 
DENVER, CO B0201 
OICY ATTN G FRE V ER 

MCDC^NELL DOUGLAS CCPP. 

5331 P n L S A AVENUE 
HUNTINGTON BEACH, CA 92647 
OICY ATTN H HE R OMAN 
OICY AT^N R HALPPIN 
OICY ATTN 0 OF AN 

MCDONNELL DOUGLAS CCPP. 

3955 LAKEWOOD BOULEVARD 
LONG BEACH, CA 90E46 
OICY ATTN M POTTER 

MERRITT CASES, INC. 

P 0 BOX 1206 
REDLANDS* CA 92373 
OICY ATTN J MERRITT 
OICY ATTN LIBRARY 

METEOROLOGY RESEARCH, INC. 

464 W WOGDRU°Y ROAD 
ALTADE M A, CA 91001 
OICY ATTN W GREEN 

MISSION PESEARCH CORP. 

P 0 DPAVF* 719 

SANTA BARBA 0 A» CA 9210? 

(ALL CLASS: ATTN: SEC OFC FOR) 
OICY ATTN C LONGMIPE 
OICY ATTN G MCCAPTCR 
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p AC IF IC-SIERRA RESEARCH CORR. 

•'1456 CL OVER Finn BLVC 
SANTA MONICA, CA 904C4 
0ICY ATTN H 3 D 00E 

PACIFTC-STERPA RESEARCH COPP. 

WASHINGTON OPERATIONS 
1401 WILSON 8L V D 
SUITE 1100 

ARLINGTON* VA 22209 
01C Y ATTN 0 GORMLEY 

PACIFICA TECHNOLOGY 
P 0 BOX 148 
DEL MAR, CA 92014 
0 IC Y ATTN R BJCRK 
01CY ATTN G KENT 
01CY ATTN TE r F LIBRARY 

PAT FL FNTFR PR IS OS s INC. 

P 0 BOX 3531 
HUNTSVILLE, AL 35R1C 
01CY ATTN M PATEL 

PHYSICS INTERNATIONAL CO. 

2700 NERCEO STREET 
SAN LEANDRO, CA 94577 
01CY ATTN L PEHRMANN 
0IC Y ATTN TECHNICAL LIBRARY 
OICY ATTN E MCORE 
01 CY ATT\' J T FC V S E K! 

OICY A T TN F SAUER 

p C 0 ASSOCIATES 
r> n n OX R655 

YAP I N A CEL PRY, CA 50291 
OICY Atjn P POP t 
OICY ATTN A KUHl 
OICY ATTN J L p T S 
01C-Y ATTN W V.' RIGHT 
OICY ATTN J CA op ENTEP 

OICY ATTN TECHNICAL INFORMATION CENTER 


RAND COPP. 

1700 mapj STREET 
SANTA MONICA, CA C 04C6 
OICY ATTN C MOW 

SCIENCE ARPLIC'T IOVS, INC 
RADIATION INSTRUMENTATION 01V 
4615 HAWKINS, NE 
ALBUO'JFROUE, NY 9?1C9 
OICY ATTN J DISHCN 
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SCIFNCE APPLICATIONS, INC. 
p 0 BOX 2351 
Li JOLLA, CA 5203<? 

0 If. Y ATTN » WILSCN 

OICY ATTN TECHNICAL LIBRARY 

OICY ATTN R SCHLAUG 

SCIENCE A?°LICAT IONS, INC. 

101 CONTINENTAL PL VO 
FL SEGUNDO, CA OC2A5 
OICY ATTN D HOVE 

SCIENCE A OD LICA T IONS* INC. 

2450 WASHINGTON AVENUE 
SAN LEANDRO, CA 94577 
OICY ATTN D BERNSTEIN 
OICY ATTN D MAXWELL 

SCIENCE APPLICATIONS, INC. 

P 0 BOX 1303 
MCLEAN, VA ??102 

01C-Y ATTN J COCKAYNE 
OICY ATTN B CHAMBERS III 
OICY A T TN M KNASEL 
OICY ATTN W LAYSCN 
OICY ATTN R SIEVERS 

SOUTHWEST RESEAPCH INSTITUTE 
P C DRAWER 2S510 
SAN ANTONIO, TX 782P4 
OICY ATTN A WENZEL 
OICY ATTN W BAKER 

SRI INTERNATIONAL 
333 PAVENSKOED AVENUE 
MENLO PARK, CA 94025 

OICY ATTN G APRAHA w SCN 
OICY ATTN LIBRARY 
OICY ATTN J COLTCN 

SYSTEMS, SCIENCE fi SOFTWARE INC 
P 0 BOX 3243 
ALBUQUERQUE mm P7198 
OICY ATTN C NEEDHAM 
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SYSTEMS, SCIENCE C SOFTWARE, INC. 
P 0 BOX 1620 
LA JOLLA, CA 92038 
01CY ATTN J BARTFEL 
0 1C Y ATTN T R INF Y 
0 1CY ATTN D C-RTNE 
01CY ATTN LIBRARY 
01CY ATTN C HASTING 
01CY ATTN K PYATT 
OICY ATTN C DISMUKES 
01CY ATTN T CHERRY 

SYSTEMS, SCIFNCE £ SOFTWARE , INC. 
1L U 00 SUNRISE VALLEY DRIVE 
RESTCN, VA 22091 
0 1CY A T TN J MURPHY 


TELEDYNE BPOUN ENGINEERING 
CUMMTNGS RESEARCH PARK 
HUNTSVILLE, AL 3S807 

OICY ATTN J PAVENSCRAFT 
OICY ATTN J MCSUAIN 


TEP°A TEK, INC. 

420 WAKAP.A WAY 
SALT LAKE CITY, LT F41C8 
01CY ATTN A ABO'I-SAYFC 
OICY ATTN LIBRARY 
OICY ATTN A JCNES 
OICY ATTN S GREEN 

TETRA TFCH, INC. 

630 N ROSEMEAD BLVD 
PASADENA, CA 91107 
01 C’Y ATTN L HWANG 

TRW DEFENSE t S^ACF SYS GROUP 
ONE SPACE PARK 
REDONDO 8EACH, CA oC?79 
OICY ATTN N LIPNER 

OICY ATTN TECHNICAL INEOPNATION CENTER 
OICY ATTN T MAZZCLA 

TRW DEFENSF £ SPACE SYS GROUP 
R 0 BOX 1310 

SAN BERNARDINO, CA 9240? 

OICY ATTN G HULCHEP 
OICY ATTN P DAI 
OICY ATTN E WCNG 


30 




UMIVFPSAL ANALYTICS, jv C . 
7740 W MANCHESTER RLVC 
° L A Y A DEL REY, C& 9C?E1 
01CY ATTN E P IEL9 


WEIMINCER ASSOC., ccnsulting engineers 

110 E 59TH STREET 
NEW YORK, NY 100 72 
01C.Y ATTN I SANOLEP 
0 1CY ATJN M PARON 


WEIDLINGER ASSOC., CCNSULTING ENGINEERS 
3000 SAND HILL PCAD 
MENLO PARK, CA 940?5 
01CY ATTN J ISENREPG 
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